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EXECUTIVE  SUMMARY 


OBJECTIVE 

Determine  the  effectiveness  of  matched-field  processing  for  source  localization  in  a  shallow- 
water  environment. 

APPROACH 

Narrowband  and  multitone  matched-field  processing  were  performed  on  vertical  line  array 
data  collected  during  the  SWellEX-1  experiment,  which  was  conducted  off  the  coast  of  San 
Diego  in  August  1993. 

RESULTS 

While  narrowband  matched-field  processing  tended  to  produce  ambiguous  source  localiza¬ 
tions  due  to  inexact  propagation  modeling,  it  was  possible  to  accurately  track  the  source  by 
simultaneously  processing  the  transmitted  tones. 

CONCLUSIONS 

Multitone  processing  was  effective  in  reducing  the  sensitivity  to  mismatch  that  can  degrade 
the  performance  of  narrowband  matched-field  processing.  Multitone  and  broadband  matched- 
field  processing  should  be  effective  tools  for  source  localization  and  tracking. 
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INTRODUCTION 


In  August  1993,  Shallow  Water  Evaluation  Cell  Experiment  No.  1  (SWellEX-1),  a  study  of 
propagation,  matched-field  processing  (MFP)  and  ambient  noise  was  conducted  in  a  shallow-wa¬ 
ter  area  approximately  5  miles  west  of  San  Diego,  CA.  Figure  1  presents  an  overview  of  the  site. 
Data  were  recorded  on  a  48-element,  88.125-m  aperture,  vertical  line  array  suspended  from  the 
Scripps  Institution  of  Oceanography’s  (SIO’s)  research  vessel  FLIP,  which  was  moored  at 
site  “O”  in  figure  1.  Acoustic  element  localization  performed  by  SIO’s  Marine  Physical  Labora¬ 
tory  put  the  bottom  element  approximately  11  m  off  the  bottom.  Water  depth  at  the  array  site  was 
approximately  200  meters,  so  the  array  essentially  spanned  the  bottom  half  of  the  water  column. 
A  source  emitting  tonals  from  70  to  755  Hz  was  towed  at  depths  from  20  to  80  m  along  an 
approximately  range-independent  path,  G  to  I,  and  along  a  strongly  range-dependent  path,  A 
to  E.  Bathymetric  measurements  and  conductivity-temperature-depth  (CTD)  casts  were  made  to 
support  the  analysis  of  the  experimental  data. 

This  report  presents  the  results  of  matched-field  processing  in  the  SWellEX-1  experiment. 
The  second  section  briefly  discusses  MFP  itself  and  presents  applications  to  data  for  the  G  to  I 
track.  The  third  section  is  a  similar  discussion  for  the  A  to  E  track.  The  last  section  summarizes 
the  findings. 


MATCHED-FIELD  PROCESSING  ALONG  G  TO  I  TRACK 


One  of  the  objectives  of  SWellEX-1  was  to  study  the  effectiveness  of  matched-field  proces¬ 
sing  as  a  means  of  source  localization  in  a  shallow-water  environment.  Originally  proposed  by 
Bucker1  in  1976,  MFP  has  been  an  active  area  of  research  in  underwater  acoustics  in  recent 
years.  Basically,  MFP  localizes  a  source  by  correlating  the  measured  sound  field  at  a  hydrophone 
array  with  that  predicted  by  a  propagation  model  for  a  series  of  candidate  source  depths  and 
ranges.  A  correlation  peak  indicates  a  high  probability  that  the  assumed  source  location  corre¬ 
sponds  to  the  true  position.  The  correlation  is  defined  as 


C(f,r,z ) 


i  j _ 

Xi?«^2i‘s^r,2)i2 

i  i 


(i) 


Here, /is  the  source  frequency,  and  r  and  z  are,  respectively,  the  range  and  depth  of  the 
assumed  source  location.  The  indices  i  and  j  refer  to  hydrophone  numbers  and  the  sums  are  over 
the  elements  of  the  array.  Si(f,rjz)  is  the  theoretical  sound  pressure  at  phone  i  produced  by  a 
source  of  frequency /located  at  range  r  and  depth  z.  This  was  calculated  using  Collins’  Finite 
Element  Parabolic  Equation  program.2  If  pi(f)  is  the  measured  sound  field  at  phone  i  and  fre¬ 
quency/,  then  the  covariance  matrix,  Rij(f),  is  calculated  as  <Pi(f)Pj*(f)>,  where  the  brackets 
indicate  a  time  average  and  the  *  indicates  complex  conjugate.  A  perfect  match  between  the  data 
and  prediction  would  yield  a  correlation  of  unity. 


1H.  P.  Bucker.  1976.  “Use  of  Calculated  Sound  Fields  and  Matched-Field  Detection  to  Locate  Sound 
Sources  in  Shallow  Water,”  J.  Acoust.  Soc.  Am.  59,  368-373. 

2M.  D.  Collins  and  E.  K.  Westwood.  1991.  “A  Higher-Order  Energy-Conserving  Parabolic  Equation  for 
Range-Dependent  Ocean  Depth,  Sound  Speed,  and  Density,”  J.  Acoust.  Soc.  Am.  89, 1068-1075. 
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The  G  to  I  track  was  a  nearly  range-independent  propagation  path.  Figure  2  illustrates  the 
bathymetry  and  sound  speed  profile  used  in  the  calculations.  The  sound  speed  profile  is  from  a 
CTD  cast  taken  at  site  G  at  the  same  time  of  day  as  the  data. 

The  data  set  analyzed  covered  the  time  period  from  1615  to  1736  on  17  August  1993.  The 
source  moved  outward  from  a  range  of  0.6  km  to  8.2  km  at  a  speed  of  approximately  3  knots. 
The  tow  depth  varied  somewhat  but  was  generally  between  80  and  85  m.  The  source  track  and 
the  position  of  the  receiving  array  in  the  water  column  are  shown  in  figure  2. 

Figure  3  is  a  single  phone  lofargram  for  the  0  to  250-Hz  frequency  band.  The  transmitted 
tones  at  70,  95,  145,  and  195  Hz  are  visible,  along  with  numerous  other  spectral  lines.  The  prom¬ 
inent  line  at  175  Hz  is  from  FLIP.  The  broadband  noise  is  due  to  a  passing  ship,  which  had  clos¬ 
est  point  of  approach  (CPA)  at  approximately  1640.  The  95-Hz  tone  is  nearly  obliterated  by  the 
noise  at  CPA.  The  signal-to-noise  ratio  was  approximately  20  dB  at  the  start  of  the  track  and 
dropped  to  5  to  10  dB  at  the  far  end. 

As  presented  above,  MFP  is  a  narrowband  technique.  That  is,  both  the  measured  and  pre¬ 
dicted  hydrophone  sound  fields  are  at  one  frequency.  A  problem  encountered  in  narrowband 
MFP  is  that  while  there  is  high  correlation  at  the  correct  source  location,  mismatch  between  the 
modeled  and  true  propagation  can  result  in  correlation  peaks  of  the  same  magnitude  at  other 
positions.  These  sidelobes  make  unambiguous  source  localization  difficult.  However,  if  the 
source  emits  several  frequencies,  Porter3  has  shown  that  simultaneous  processing  of  the  tones 
results  in  improved  performance.  Essentially,  each  frequency  has  a  correlation  peak  at  the  true 
source  location,  but  the  positions  of  the  sidelobes  vary  with  frequency,  so  combining  the  fre¬ 
quencies  will  lower  the  sidelobe  level  relative  to  the  true  peak.  While  there  are  numerous  ways 
in  which  the  frequencies  can  be  combined,  a  simple  incoherent  averaging  scheme  was  found  to 
be  effective.  For  each  assumed  source  location,  the  correlation  in  equation  1  is  calculated  sepa¬ 
rately  for  each  of  the  M  tones,  and  these  correlations  are  then  averaged  to  yield  the  final  result: 

C(r’z)  =  A/XC(^’r’z)  (2) 

/ 

To  illustrate,  the  top  four  frames  in  figure  4  are  the  matched-field  range-depth  surfaces  for  70, 

95, 145,  and  195  Hz,  respectively,  at  a  time  when  the  true  range  was  5.9  km  and  the  depth  was 
85  m.  The  multiple  correlation  peaks  in  these  plots  make  it  difficult  to  say  with  much  confidence 
where  the  source  actually  is.  However,  if  these  narrowband  results  are  averaged  as  in  equation  2 
with  M  =  4,  the  bottom  plot  is  obtained.  The  correlation  maximum  at  a  5.9-km  range  and  a  84-m 
depth  agrees  well  with  the  true  source  location,  and  the  localization  is  fairly  unambiguous,  in 
that  the  sidelobes  are  considerably  weaker  than  the  main  peak. 

To  see  how  well  the  source  could  be  tracked,  the  81 -minute  data  set  was  treated  as  a  set  of 
1-minute  segments,  overlapped  by  50  percent.  For  each  of  these  segments,  narrowband  and 
frequency-averaged  surfaces  like  those  in  figure  4  were  generated,  and  for  each  surface,  the 
range  and  depth  at  which  the  maximum  correlation  occurred  was  saved.  Figure  5  summarizes  the 
results.  The  left  two  frames  plot  the  range  and  depth  of  the  peak  correlation  vs  time  for  each  of 


3M.  B.  Porter.  1994.  “Source  Tracking  in  the  Hudson  Canyon  Experiment,”  J.  Acoust.  Soc.  Am.  95,  2981. 
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frequency-averaged  surfaces  like  those  in  figure  4  were  generated,  and  for  each  surface,  the 
range  and  depth  at  which  the  maximum  correlation  occurred  was  saved.  Figure  5  summarizes  the 
results.  The  left  two  frames  plot  the  range  and  depth  of  the  peak  correlation  vs  time  for  each  of 
the  individual  frequencies.  There  are  a  considerable  number  of  false  localizations  where  inexact 
modeling  has  resulted  in  a  sidelobe  having  higher  correlation  than  that  at  the  true  source  loca¬ 
tion.  The  right  two  frames  in  figure  5  are  similar  plots  for  the  frequency-averaged  case.  For  com¬ 
parison,  the  solid  lines  are  the  true  range  and  depth  vs  time.  Frequency-averaging  has  been  effec¬ 
tive  in  reducing  the  sidelobes  and  consequently  the  source  depth  and  range  tracks  are  accurately 
reproduced. 

Rather  than  plot  the  matched-field  output  as  a  function  of  range  and  depth  at  a  fixed  time  as 
in  figure  4,  we  can  fix  the  source  depth  and  plot  the  output  as  a  function  of  range  and  time.  Fig¬ 
ure  6  presents  narrowband  and  frequency-averaged  range-time  plots  for  a  source  depth  equal  to 
the  average  tow  depth,  82  m.  The  source  track  is  clearest  in  the  frequency-averaged  plot, 
although  the  145-Hz  and  195-Hz  narrowband  results  are  also  fairly  unambiguous.  As  noted  ear¬ 
lier,  the  95-Hz  tone  was  degraded  by  the  broadband  source.  The  70-Hz  tone  has  the  highest  cor¬ 
relation,  but  also  has  the  highest  sidelobes. 

Two  other  data  sets  from  the  G  to  I  track  were  analyzed  in  the  same  manner  as  above.  Only 
the  final  results  will  be  presented.  From  1745  to  1910  on  17  August,  the  source  was  towed 
inward  from  I  to  G  at  a  depth  of  approximately  80  m.  Figure  7  is  a  plot  of  the  range  and  depth  of 
the  correlation  peak  vs  time  for  this  data  set.  Finally,  from  2205  to  2340  on  17  August,  the  source 
was  again  towed  from  I  to  G,  but  this  time  at  a  shallow  depth  of  20  m.  Figure  8  summarizes  the 
narrowband  and  frequency-averaged  results.  In  these  two  data  sets,  as  in  figure  5,  frequency¬ 
averaging  has  had  the  desired  effect  of  suppressing  the  sidelobes  so  that  reliable  source  tracking 
is  possible. 

MATCHED-FIELD  PROCESSING  ALONG  A  TO  E  TRACK 

In  contrast  to  the  relatively  benign  G  to  I  track,  the  A  to  E  track  was  strongly  range-depen¬ 
dent.  Figure  9  illustrates  the  propagation  environment.  The  sound  speed  profiles  are  from  CTD 
casts  taken  at  the  indicated  locations. 

From  0956  to  1129  on  17  August,  the  source  was  towed  outward  at  an  average  depth  of  21  m 
from  an  initial  range  of  0.6  km  to  a  final  range  of  9.2  km. 

Figure  10  is  a  lofargram  for  this  data  set.  The  transmitted  tones  at  70,  95, 145,  and  195  Hz 
are  identified.  The  signal-to-noise  ratio  ranged  from  15  to  20  dB  at  point  A  to  approximately 
5  dB  at  point  E. 

Figure  11  presents  narrowband  and  frequency-averaged  range-depth  surfaces  at  a  time  when 
the  range  was  6.2  km  and  the  depth  was  21  m.  The  narrowband  results  are  very  ambiguous.  The 
correlation  peak  for  the  frequency-averaged  case  is  at  a  6.2-km  range  and  a  24-m  depth,  in 
agreement  with  the  actual  location,  but  the  sidelobes  are  only  slightly  lower  than  the  main  peak. 
In  general,  frequency-averaging  was  not  as  successful  in  suppressing  the  sidelobes  on  this  track 
as  it  was  on  the  G  to  I  track. 
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Figure  12,  generated  in  the  same  manner  as  figure  5,  plots  the  range  and  depth  of  the  peak 
correlation  vs  time.  The  solid  lines  are  the  true  range  and  depth  tracks.  Although  the  frequency- 
averaged  plot  shows  some  false  localizations,  the  track  is  still  fairly  well  defined.  It  is  a  marked 
improvement  over  the  narrowband  results,  where  there  is  hardly  a  hint  of  the  source  track. 

Finally,  figure  13  is  a  set  of  narrowband  and  frequency-averaged  range-time  plots  for  a 
source  depth  equal  to  the  average  tow  depth,  21  m.  The  source  track  is  better  defined  in  the 
frequency-averaged  plot  than  in  any  of  the  narrowband  plots,  although  none  of  the  tracks  are  as 
clear  as  their  counterparts  in  figure  6  for  the  G  to  I  track. 

CONCLUSIONS 

The  idea  of  matched-field  processing,  localizing  a  source  by  comparing  the  measured  and 
modeled  propagation,  is  appealing.  In  practice,  however,  imprecise  knowledge  of  environmental 
parameters  can  seriously  degrade  MFP’s  performance.  Through  analysis  of  data  collected  on  the 
SWellEX-1  experiment,  it  was  found  that  frequency-averaging  is  an  effective  means  for  reducing 
the  sensitivity  to  modeling  uncertainty,  and  makes  unambiguous  source  tracking  possible.  While 
the  scenario  studied  here,  four  widely  spaced  tones,  is  unlikely  to  occur  in  practice,  the  results 
indicate  that  MFP  would  be  effective  for  true  broadband  sources. 
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Figure  1.  SWellEX-1  test  area,  August  1993. 
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Figure  2.  Track  G  to  I — bathymetry,  sediment  thickness,  and  sound  speed  profile. 
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Figure  3.  Track  G  to  I,  lofargram  for  Channel  24. 
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Figure  4.  Track  G  to  1, 1710  to  1711,  true  range  =  5.9  km,  true  depth  =  85  m. 


9 


NARROWBAND  RANGE  TRACKS  FREQUENCY-AVERAGED  RANGE  TRACK 
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Figure  6.  Track  G  to  1, 1615  to  1735,  source  depth  =  82  m. 
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Figure  9.  Track  A  to  E — bathymetry,  sediment  thickness,  and  sound  speed  profiles. 
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Figure  10.  Track  A  to  E,  lofargram  for  Channel  24. 
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Figure  11.  Track  A  to  E.  1056  to  1057.  true  range  =  6.2  km.  true  depth  =  21  m. 
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Figure  13.  Track  A  to  E,  0956  to  1126.  source  depth  =  21  m. 
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